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Abstract 

Hie angles of propagation of the wave fronts 
associated with duct modes are derived for a cylin- 
drical duet with soft walls (acoustic suppressors) 
and a uniform ateady flow. The angle of propaga- 
tion with rcapcct to the radial coordinate (angle 
of Incidence on the wall) is shown to bo a better 
correlating parameter for Che optimum wall Imped- 
ance of spinning modes than the previously uoed 
mode cutoff ratio. Doth the angle of incidence 
upon the duct wall and the propagation angle with 
respect to the duct axiB are required to describe 
the attenuation of a propagating mode. Using the 
modal propagation angles, a geometric acoustics ap- 
proach to suppressor acoustic performance was de- 
veloped, Results from this approximate method were 
compared to exact modal propagation calculations to 
check the accuracy of the approximate method. The 
results arc favorable except in the Immediate vicin- 
ity of the modal optimum impedance where the approx- 
imate method yields about one-half of the exact 
maximum attenuation. 


Introduction 

The angles of propagation for the wave fronts 
making up on acoustic mode arc derived for ducts 
with soft wslla (acoustic liner). These are the 
anglea which the velocity vector normal to the 
local wave front makes with tile three coordinate 
axes. The angles associated with the group veloc- 
ity vector are also derived. Tills effort is an ex- 
tension of the work of Ref, 1 in which both approx- 
imate and exact wave propagation angles were de- 
rived for a cylindrical duct with uniform steady 
flow. In Ref, l, only the hardwall results were 
given and the emphasis was on the relationship of 
the propagation anglea in Che duct and the far- 
ficld radiation pattern. This paper will consider 
ducts with soft walla and Che relationship of the 
propagation angles to acoustic liner performance. 

In Refs. 2 to 4 Che optimum acoustic impedance 
and the maximum possible attenuation for an acous- 
tic mode were shown to be intimately connected to 
and correlated by the mode cut-off ratio. The mode 
cut-off ratio Is a somewhat abstract parameter, 
particularly for soft wall ducts, and an explana- 
tion of its effectiveness as a correlation parame- 
ter was that it was connected to the angle of in- 
cidence upon the liner wall. This paper formalizes 
this connection and shows that angle of incidence 
is the more fundamental parameter and in fact cor- 
relates Che optimum impedance values better than 
does the cut-off ratio. 

The maximum possible attenuation for a mode 
was also shown to be approximately correlated by 
the cut-off ratio. However, the collapse of these 
attenuation curves was not as good as desired and 
some large errors could occur particularly near 
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mode cut-off when using this correlation. A sim- 
ple geometric approach nimllar to that of Ref, 5 
was tried. Hits approach considers that the num- 
ber of encounters that the wave has with the soft 
wall is determined by the wave angle with respect 
to the duct axis end by the duct size. The acous- 
tic power remaining after each encounter with Che 
wall is rolatod ta angle of Incidence and the wall 
impedance properties, When these two concepts are 
combined in a simple model, the acoustic power at- 
tenuation for the mode can bo calculated. The re- 
sulting equation is shown to be exact for a rect- 
angular duct without Mach number and for only 
slightly soft walls. This approximate method is 
extended to Include cylindrical ducts and unifotm 
duct flow. Results from the approximate model are 
compared with the exact attenuation results for 
several conditions Including acoustic impedances up 
to the optimum for both rectangular and cylindrical 
ducts with flow. 

The purpose of this paper is to provide the 
foundation for a simple geometric acoustics ap- 
proach which may he useful for complex suppressor 
geometries end multimodal noise sources for v.'h < <'h 
exact methods are inefficient and extremely diffi- 
cult. 


Symbols 

c speed of sound, or the magnitude of the 

wave velocity vector normal to a plane 
wave front, m/sec 

C£ components of vector e, i ° x, r, C>, m/sec 
cr resultant velocity vector, ■'« 

D duct diameter, m 

f frequency, Hz 

H rectangular duct height, m 

Ip axiol acoustic intensity leaving end of 

suppressor, N/m sec 

Ip axial acoustic intensity entering suppressor, 
N/m sec 

I x axiol acoustic intensity of reflected wave 

. off soft wall, N/m sec 

Ix axial acoustic intensity of incidence wave 

on soft wall, N/m sec 
l 

J m Bessel function of the first kind of 
order m 

j rectangular duct mode number 

k wave number, to/c, m"^ 

k r radial wave number, m - ^ 

kp], combined radial-circumferential wave num- 

ber (a/r 0 ), m' 1 
k x axial wave number, m"* 

kj, circumferential wave number, m"^ 

L suppressor length, ra 

M Mach number of uniform steady duct flow 

M 0 Mach number of uniform outside of boundary 

layer 

m -pinning mode lobe number 

P _cou8tic pressure, N/m^ 

t radial coordinate, m 

r 0 duct radius, m 



R magnitude of complex eigenvalue 

t time, see 

x axial coordinate, m 

Ax axial distance between encounters of acous- 

tic wave with duct walls, m 
y transverse rectangular coordinate normal to 

duct walls, m 

z transvorso rectangular coordinate parallel 

to duct walls, m 
a complex eigenvalue, Re^T 

p real part of complex quantity, ace Eqs, (12), 

(13) 

6 boundary layer thickness for l/7th power 

velocity profile, m 

e dimensionless boundary layer thickness, 6/r 0 

T) frequency parameter, fD/c or fli/C 

0 specific acoustic resistance 
optimum specific acoustic resistance 

p radial mode number 

£ mode cut-off ratio, see Eq. (22) 

a attenuation or damping coefficient, see 

Eq. (7) 

T propagation coefficients, see Eqs. (7), (8) 

<t> circumferential coordinate, radians 

9 phase of eigenvalue a, radians 

<p r wave front propagation angle relative to 

radial coordinate, angle of incidence on 
the wall, deg 

q> x wave front propagation angle relative to 
axial coordinate, deg 

tPy wave front propagation angle measured from 
y coordinate, deg 

qfy wave front propagation angle relative to 

circumferential coordinate, deg 
X specific acoustic rcsctsncc 

Xn, optimum specific acoustic reactance 

resultant axial propagation angle in duct, 
deg 

resultant propagation angle measured from 
radial coordinate, dm? 

V,j, resultant propagation angle measured from 

circumferential coordinate, deg 
(0 circular frequency, radlans/sec 

Subscripts : 

1 equal to X, r, or $ 

op mode indices, m lobe p radial order 
R indicates real port of complex quantity 

r 0 evaluated at outerwall 


Development of the Propagation Angies 


■ The initial development of the expressions for 

■ the propagation angles can proceed exactly ns in 
Ref, 1. The wave equation in a circular duct with 
a steady flow can be expressed as 


I (1 - H 2 ) — - 
ox 2 


2M a 2 P c£p 

c dxdt Sr z 


, i 2 p , 
+ x 5 ? + 


_L s^p „ _i_ a^p 

r 2 &J> 2 c 2 3t 2 


a> 


(where x, r, and <J> are the usual cylindrical 
coordinates, t is time, M is Mach number, c the 
.speed of sound, and p is the acoustic pressure. 
The solution to Eq. (1) for the pressure is, 



where m is the number of lobes for the spinning 
mode (the transverse order), p is the radial order 
number of the mode, kx,ny is the axial wave num- 
ber, ai is the circular frequency, r Q is the outer 
wall radius, is the eigenvalue of the iqt 

mode, and J m is the Bessel function of the first 
kind and order m. 


Wave Humbcrs and the Phase Velocity 
Propagation Angles 

For brevity, the iqi subscripts will be de- 
leted from hero on, and it will be understood that 
s single (but quite general) mode is being con- 
sidered. When Eq, (2) is inserted into Eq. (1) 
the result is, 

2 

(r^) + k x H k2 * 2l,lkk x + M 2 k 2 » (k - Mk x ) 2 (3) 

The first term in Eq. (3) is a combined radial- 
transverse wave number squared which could be de- 
noted by, 


Tills combination of wave 'lumbers is the causi? of 
the trouble in defining a ime of the propagation 
angles in cylindrical ducts, 


The i fc li propagation angle can be shown (for 
hurdwall ducts), 1 to be given by, 


cos (p£ 


ki 


k l 


4 -t) 


(5) 


where i = r, $, or x, and (p^ is the angle be- 
tween the normal to the wave front and the i tk 
coordinate axis, and the denominator of Eq. (S) 
was obtained from Eq. (3). The angles cp x , (p r 
and cp^> can be seen in Fig. 1. The vector c is 
normal to the local wavefront. The other angle 
set and the resultant vector cp due to addi- 

tional convective effects nlso shown in Fig. 1 will 
be discussed later. 


Before proceding further, one must recognize 
that, due to the soft wall boundary condition, the 
quantities a and kj^ in Eq. (3) are complex num- 
bers Only the real parts of the wave number are 
related to propagation wave numbers for use in 
Eq. (5) while the imaginary parts represent damp- 
ing terms. The axial wave number can be rewritten 
as, 


U x - k(x - in) 

and this alon». with Eq. (3) yields, 


-iM + i 


o + iT = 


j/i - a - « ; )(^) Z 


K* 


( 6 ) 


(7) 


and thus , 


2 


( 8 ) 


-M + 


„ 1/2 

r 2 [ l * <'■ - m2 >© 2 cqb 2 * + / l ^ - M2) (^) r * 2(1 ' 


1 - m2 


where R and q> are the amplitude and phase of 
the complex eigenvalue a and tho frequency param- 
eter is, 





° ire 


ft) 

c 


O) 


Equation (8) is obtained by solving for the imagi- 
nary part of Eq. (7) . 

To obtain the angle of incidence ac the linor 
wall (<p r ) in a simple form, the assumption is made 
that the combined radial- transverse complex wave 
number can be split as follows, 



Group Velocity Propagation Angles 

Equations (15) to (17) describe the orienta- 
tion of the wove front at the outer wall of the 
circular duct and can be seen to include the duct 
Mach number through tho propagation coefficient T, 
There is an additional Mach number effect os shown 
in Fig. 1 which causes a drift in the axial veloc- 
ity and produces a resultant vector (eg) of propa- 
gation which is not normal to the wavefront. This 
vector is associated with the group velocity*- and 
the angles of propagation associated with eg will 
be of use in a later discussion. These angles can 
be derived from the geometry shown in Fig. 1 using 
tho face that tho vector cM alters only thu x 
component of velocity but not Cf or c$ . These 
angles can be expressed ns, 


Tills assumption was used in Ref, 1 for hordwnll 
dusts and was found to provide adequate results 
near the outer wall when compared to on exact solu- 
tion using Hankcl functions. The real part of 
k r (k r( g) for use in Eq. (5) can then be obtained as, 



M + cos <p x 

cos “ — 1 (18) 

yl + M 2 + 2M cos {p x 

cos qj_ 

cos iff,. « . -= (19) 

V 1 + M 2 + 2M cos qi x 

and 


( 11 ) 

Tho propagation angles will be calculated only at 
the outer wall (r = r 0 ) whore the approximations 
used here ore reasonably accurate.*- Let, 

kr - R ) t 0 - S (12> 


thus 


P 


— (r 2 cos 2<p - m 2 + - 2m 2 R 2 cos 2ip + in*) 

vS 


(13) 

The circumferential wave number (at the outer woi.) 
as seen from Eq. (10) is, 


kq, 


jn mk 

r D ° in) 


(U) 


which is real unlike the kx and k r . 

The real parts of Che wave numbers can now be used 
in Eq, (5) to obtain, 


cos <p x 


cos tp r 


+ (n) 2 * (S) 2 


(15) 


(16) 


cos 


irqyv + (^) 2 + (^’ 




(17) 


cos 5|fj, 


cos 

VT + M 2 + 2M cos q> x 


( 20 ) 


Correlation of Optimum Impedance 
with Incidence Angle 

Hie angle of incidence on the soft wall 
(Eq. (16)) was used to correlate some of the opti- 
mum impedance calculations reported in Ref. 4. 

These calculations were made for a cylindrical 
duct with a l/7th power boundary layer and a uni- 
form core flow outside of the boundary layer. 
Standard modal solutions wore used within tho core 
flow and were coupled to the wall impedance by a 
Runge-Kutta integration through the boundary layer. 
These correlations are shown in Figs. 2 and 3. 

These calculations include several lobe numbers 
(m) and the entire range of radial modes from well 
propagating to near cut-off. Comparison of FigB. 2 
and 3 with Figs, 12 and 13 of Ref. 4 show that the 
angle of incidence correlates the modal optimum 
impedance somewhat better than the mode cut-off 
ratio did in Ref. 4. Using cut-off ratio^ the 
first radial mode for m “ 5 showed some deviation 
from the correlation while m =■ 9 and 16 had quite 
pronounced deviations. This deviation does not 
appeor when incidence angle is used. This improve- 
ment can be most easily explained by using the 
hardwall incidence angle given in Ref. 1 which was 



vnzi 

&(i - m y/\ - i/s 2 ) 


where the cut-off ratio £ is given by 


( 21 ) 
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kr n 


a\/l * H z 


( 22 ) 


For higher radio! mod es (a. » bi) or for small m 
the terra VI (m/a)2 approaches unity, Equa- 
tlon (21) then shows that for constant Mach number 
M (as used In Figs. 2 and 3) the angle of Incidence 
(<pj:) Is o function only of cut-off ratio and cither 
parameter Is sufficient to correlate the optimum 
impedance. However, for large m and low radial 
order a » m, the term Vl - (m/a,)2 in Eq. r 21) In- 
fluences the incidence angle and causes Scatter In 
the correlation as shown in Ref. 4. At sufficient- 
ly high lobe numbers Che splitting of the radial 
and transverse wave numbers as used In Eq. (10), a 
rectangular coordinate approximation to the cylin- 
drical coordinate system, will Iobg validity and 
some scatter in the corrleatlon will again be in- 
troduced. This latter problem can probably be 
overcome by extending the exact Hankcl function 
solutlona of Ref. 1 to soft walls. 


The optlmu' impedance correlation shown in 
Figs, 2 and 3 Wu.'e for propagation calculations 
containing a boundary layer. The angle of Inci- 
dence calculated in each case wsb valid only in the 
uniform flow region outside the boundary layer. 

A complete geometric approach to duct propagation 
must contain the refractive effects on the acoustic 
wave as it passes through this boundary layer. 

This is emphasized by the calculated optimum imped- 
ances shown in FigB. 4 and 5 for three boundary 
layer thicknesses. Note that for onglea of inci- 
dence (cpf) above 65° there is a dramatic effect of 
boundary layer thickness upon optimum resistance 
(Fig, 4) while the reactance (Fig. 5) is affected 
down to about 50°. For a given set of propagation 
angles outside of the boundary layer there is a 
unique optimum impedance but the boundary layer 
acts as an impedance transformer between the uni- 
form flow region and the wall. Fortunately for 
modes neorei. to cut-off, the wave motion is not 
directed into the velocity gradients (mainly radial 
or circumferential) and the boundary layer has much 
less effect upon the propagation of these modes. 

An alternate correlation of one set of calcu- 
lations (e ° 0.005) in Fig. 5 is shown in Fig. 6. 
Here the resultant incidence angle (^ r ) is used 
rather than tp r - Recall that V r contains the 
axial velocity drift effect and is associated with 
the group velocity. Some of the scatter in the 
points nearer cut-off (small incidence angle) may 
have been reduced by using and the numerical 

value of incidence angle indicates more of a normal 
incidence than does (p^C^r <p r )* The fi seem to 
have more physical significance than the cpj.. For 
example in Ref. 1 it wa3 shown that # x - 90° at 
mode cut-off while cp x does not. 

From the above discussion it is thus proposed 
that the angle of incidence of tho wove fronts on 
| the lined wall is a more fundamental and obviously 
■more easily visualized parameter than the mode cut- 
off ratio for correloting liner optimum impedance. 
However, the mode cut-off ratio is still considered 
■ to be an extremely useful parameter for liner de- 
sign purposes with the previous discussion serving 
; to explain the reason for its effectiveness and 
pointing out its limitations. 


Estimation of Attenuation Using 
Propagation Angles 

Tho previous section dc«it only with the con- 
nection of optimum impedance with wove ongle of in- 
cidence, To complete the discussion some estimate 
of sound power attenuation muBt be made. An esti- 
mate of liner performance can bo obtained using the 
method of Ref. 3 along with the improved impedance 
correlations of Ref. 4. Both of these references 
are based upon mode cut-off ratio. However, in 
this section o geometric approach similar to those 
of Refs. 5 and 6 will be uacd, Tills ray acoustics 
approach using the propagation angles derived in 
the previous section clearly illustrates the impor- 
tant aspects of the Bound propagation phenomena in 
lined ducts which ore usually obscured in the phys- 
ical acoustics approach. It is all the more rotis- 
factory if usable rosultB for liner design can be 
obtained. The first example will consider tho sim- 
ple two dimensional rectangular system without flow 
to illustrate th principles and provide motivation 
for an cxpoctati/O of useful results. The analysis 
will then be extended to include uniform flow in o 
rectangular duct. Using these nrlnciplcB an acous- 
tically lined cylindrical duct will then be con- 
sidered , 


Rectangular Two Dimensional Duct. M ■ 0 


Figure 7 ohows a rectangular duct without flow 
with a transverse mode propagating ond bouncing be- 
tween tho soft walls. The angle of incidence be- 
tween the velocity vector normal to these waves and 
the y coordinate is tpy while the angle of the 
waves with respect to the duct axis is ep x , The 
axial distance between encounters between the wave 
and the soft wall is, 


Ax 


H 

tan tp x 


(23) 


where H is the duct height. The total number of 
walL bounces in a duct of length L Is then 


N 



5 tan tp x 


(24) 


The ratio of axial flux after a wall bounce com- 
pared to the incident wave is given by 

i; <e 2 + X^)cos^q) y - 20 cos <p y + i 

-f - — r ; ~ (25) 

I x ( 02 + X' i )COS / tpy + 20 cos tp y + 1 

where 0 and X are the wall specific acoustic 
resistance and reactance. Note that Eq. (25) in- 
volves only the wall properties and tho angle of 
incidence. Equation (25) was derived in the same 
way as the absorption coefficient is conventionally 
derived. A plane wave is assumed to be incident 
upon a plane absorber and the angle of reflection 
is equal to the angle of incidence. The absorption 
coefficient is 1 - IJj/ljj. 


After the N bounces in the duct the ratio of 
the final to initial axial acoustic Intensity is, 


h 



(26) 
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Since only one mode to considered here the attenua- 
tlon con be obtained from the Intensity ratio os, 

Ip I x 

AdD •> 10 Log ~ - ION log — (27) 

!t I + 

*• A x 

Now consider the often studied problem of 
nearly hard walls which Is quite a good approxima- 
tion to the exact physical acoustics solutions ex- 
cept In Che vicinity of the optimum impedance. 

With only slightly soft walls, 


ond 



46 

COS tpy(6^ + X 2 ) 



-J*0 

In 10 cos <p y (6 2 + X 2 ) 


(28) 


(29) 


Equations (28) and (29) were derived using Eq. (2"}) 
with the assumption that impedance is very large. 

To the firBt order of approximation, the axial 
propagation angle and the incidence angle for the 
jth transverse mode ore given by 


tan !p x 


J 

n Vi - (J/n ) 2 


(30) 


and 


cos q>y » 


J 

~) 


(31) 


Equations (15) and (16) were used with y replac- 
ing r, m ■* 0, and the nearly hardwoll approxima- 
tion was used to evaluotc the eigenvalues. When 
Eqs . (24), (29), (30), and (31) are used in Eq . (27) 
the result Is, 


AdB » 


-17.37 3 (L/H) 

Vi - <J/n ) 2 ( 02 + X 2 ) 


(32) 


This equation for attenuation is exactly the some 
as would be obtained from the physical acoustics 
approach to the same nearly hardwoll problem. Thus 
ot least for rectangular ducts without steady flow 
and with wall impedance not in the vicinity of the 
optimum impedance, the ray acoustics approach using 
the angles of propagation and the acoustic flux 
loss per wall bounce gives very precise ottenua- 
' tions . 


Rectangular Duct with Uniform Flow 

i 

With a finite Mach number in the duct Eq. (27) 
I remains valid but the number of encounters of the 
.wave with the wall (Eq, (24)) and the ratio of the 
1 incident to reflected axiol intensity (Eq. (25)) 
[must be modified. The number of bounces must now 
! include the axial propagation angle associated with 
: the group velocity (^x). Equation (24) thus be- 
I comes , 

I 

, N - J tan ( 33 > 

where Eq. (18) is used for ^x" 


The ratios of axial Intensities were dorlvcd 
as in the previous section with the addition that 
uniform flow must be considered in the wove equa- 
tion and the continuity of displacement raunt bo 
used at the soft wall. The derivation proceeds ex- 
actly as In Ref. 7, except thnt the relation be- 
tween the wave numbors and the propagation angle 
was incorrectly expressed In Ref. 7 (effect of H 
on angles Ignored), The axial Intensity ratio (re- 
flected to incident) is, 

“ [(6 2 + X 2 )coa 2 qiy(l + M sin <py) 2 
*x 

- 26 con q> y (l + M sin tyy) + 

[(6 2 + X 2 )cos 2 'p y (l + M sin ?y) 2 

+ 26 cos (py(l + M sin qjy) + lj (34) 

Note that Eq. (34) reduces to Eq. (25) when M * 0. 

For the approximate colculations to be shown 
later, Eqs. (15) and (15) were used with y re- 
placing r ond with m = 0. If a three dimension- 
al rectangular duct calculation is desired with a 
transversly spinning wav<i considered, m could be 
retained and the 2 coordinate would Just replace 
the $ coordinate. 



Cylindrical Duct with Uniform Flow 

For a cylindrical duct with uniform steady 
flow the results of the previous section can be 
used directly as a first approximation. Of course 
<p r must be used instead of cpy in Eq. (34). As 
shown in Ref, l, the approximations leading to 
Eqs. (15) to (17) are valid only near the outer 
wail because of curvature effects in the cylindri- 
cal duct. Equation (34) represents the physics of 
reflection and absorption which occur at the wall 
and it should thus remain valid. However, the num- 
ber of encounters of the wave with the wall as ex- 
pressed by Eq. (33) involves the propagotion pro- 
cess throughout the entire depth or radius of the 
duct. For a rectangular duct ^x does not c bsnge 
with distance from the wall and Eq. (33) is thus 
valid. For a cylindrical duct, as discussed in 
Ref, 1, is a function of radius and is alwnyB 
equal to 90° at the duct centerline. Thus using 
the value of # x ot t ' ,e wall, as done here, will 
underestimate the number of bounces somewhat, and 
the approximate calculation of attenuation wil! be 
less than the exact attenuation. Some integrated 
value of <r x muaC be used rather than the value at 
the wall, but this has not yet been accomplished. 


Comparison of Approximate and 
Exact Attenuations 


A series of exact single mode propagation col- 
culations were performed as in Ref, 8, These cal- 
culations were made for both rectangular and cir- 
cular ducts. In each case a constant damping coef- 
ficient (a) was used and the propagation coeffi- 
cient (t) was varied to produce a closed loop in 
the wall impedance plane (see Fig. 4, Ref. 8). 
Several points were used along this constant atten-r 
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uatlon contour and the values of the variables 
needed to calculate the approximate attenuation 
(0, X, -J, T, R, tp) wore recorded, 'rims at those 
points the axact and approximate attenuations could 
be compared. The damping value was then Increased, 
repeating the same procedure, until the damping 
contour collapsed to a point representing the maxi- 
mum possible attenuation for the mode used. The 
results are presented In the following sections, 


Rectangular Duct 

The ratio of the attenuation by the approxi- 
mate calculation using the geometrical acoustics 
approach to that of the exact physical acoustics 
calculations are shown in Fig, 8 for zero steady 
flow, Tiie abscissa represents the degree of atten- 
uation used for the calculations with the value of 
unity representing tho maximum possible attenuation 
for this r tde, the second symmetric mode. The sym- 
bols ropresent tho average attenuation for tho sev- 
eral points on a given exact damping contour while 
the bars represent the spread of these calculations , 
The averages are not unique, particularly where a 
large spread is present, since they depend upon the 
particular choice of points along the damping con- 
tour, Sevorol observations can be made from Fig. 8 . 
As damping is reduced (far from tho optimum imped- 
ance) the approximate and exact attenuations are 
seen to be in better agreement with very little 
scatter. This agrees with the comments mode per- 
taining to the limiting Eq, (32) which states that 
the approximate equation is exact far from the 
optimum impedance. Up to about one-half of the 
maximum attenuation the average attenuation around 
a contour is very good although considerable scat- 
ter occurB , As the optimum impedance iB approached 
the approximate calculation giveB only about one- 
half the attenuation as the exact calculation. At 
first glance, the 50 percent error possible here 
both at intermediate and high damping may appear 
unacceptable. However, the upper half of the damp- 
ing range represents a very small portion of the 
impedance plane which would be very hard to realize 
precisely in a real application. Furthermore, in a 
multimodal excitation situation, where approximate 
techniques would be most valuable, a given wall im- 
pedance would be near to only a small fraction of 
the modal optima and errors involving only o few 
modCB (assuming near cquipartion of energy) would 
not seriously affect the final result. Also the 
: apparently correct averaging where considerable 
scatter occurs would help the final multimodal re- 
sult, The reason for the scatter observed will be- 
come more clear in a later discussion. 

In Fig. 9 results are shown for the same mode 
; as in Fig. 8 except that there 1 b a steady inlet 
! flow of H “ -0.4. The results are much tho same 
' as In Fig. 8 with the law damping points showing 
j agreement between approximate and exact attenuo- 
- tlons which indicates that the Mach number effects 
j were handled properly, Several points are shown at 
I the Intermediate damping point where the large 
[ scatter occurs. Two points are seen to be about 
' 50 percent high while the bulk of the points on 
this damping contour are within about 20 percent of 
the exact attenuation. 

To study the scattering of the approximate 
attenuation calculations, s complete scan over all 
values of propagation coefficient T was made 
□long n constant damping coefficient contour of 


o = 0,02. Tho approximate to exact attenuation 
ratio is shown in Fig. 10 along with the incidence 
angle tfy and the magnitude of the eigenvalue R, 
Thlo scan Involves several closed loops In the im- 
pedance plane and of course cncomposaca several 
modes, in Fig. 10(a) the scatter in approximate 
attenuation Is scan to be s smooth oscillatory 
function, The scatter is seen to bo small for the 
hlghar modes (mode number noted near top) and 
increaacs for the lower modes, The value of damp- 
ing used for this calculation could not be attained 
by the first mode. There is no error at the bor- 
ders of each mode and the average within a mode can 
be seen to be about unity in the attenuation ratio. 
If a cross plot wore shown between Figs, 10(a) und 
(c) it would be seen that the scatter Is cyclic 
with eigenvalue amplitude. Tho scatter seems to 
be least for the higher, neor cut-off modes nearer 
to normal incidence (small q>y) and to increase 
with angle of Incidence. Tho^ generality of these 
resultB have not boon tested. 


Cylindrical Duct 

Figure 1L shows comparisons of approximate and 
exact attenuations for a cylindrical duct for the 
aecond radial of a seven lobed spinning mode. The 
scatter is seen to be extremely small when compared 
to Figs. 8 and 9, This 1 b probably due to the high 
frequency (q = 20 versus 5) au compared to the 
rectangular duct results. The same roll-off at 
high damping is Been here as in the rectangular 
duct. A difference occurs at low damping with the 
asymptote being about 0.9 rather than unity as pre- 
viously seen, Tills is no doubt duo to the under- 
estimation of the number of encounters of the wave 
with the wall due to using the axial propagation at 
the wall rather than an integrated value over the 
duct. Recoil that axial propagation angle Is con- 
stant over a rectangular duct but increases to 
90° at the center of n cylindrical duct. 1 


Concluding Remarks 

The angle of incidence has been shown to be a 
more fundamental and more illuminating parameter 
than the mode cut-off ratio for correlating modal 
optimum impedances. The liner performance has been 
shown to involve both the angle of incidence and 
the axial angle of propagation along with the Usual 
parameters of wall impedance and duct geometry. 

The use of these propagation angles in e ray ucous” 
tics propagation approach shows great promise al'* 
though some additional effort 1 b still needed, to 
perfect this approach. If the ray acoustics ap- 
proach con be perfected it would be an extremely 
powerful tool to solve problems in soft wall duct 
configurations which can not be presently handled 
except by costly numerical procedures. An example 
is refraction effects in variable area ducts with 
flow provided that the reflection effects are not. 
dominant. Another important point is that the ex- 
pression of modal properties in terms of propaga- 
tion angles is nearly equivalent to an expression 
in terms of cut-off ratio. The same statements can 
thus be made in terms of propagation angles as were 
mode in Refs. 2 to 4 in terms of cut-off ratio. 
Modes with nearly equivalent propagation angles 
will behave similarly In an acoustic liner. Agroup 
of modes with similar propagation angles can thus 
be treated as an entity without the need to distin- 
guish between the individual modes. Tills is espe- 
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dally significant if a large number of propagating 
modes arc available and carrying significant acous- 
tic power. If a direct measurement of acoustic 
power In the duct can be made as a function of In- 
crements in the propagation angles, then a signifi- 
cant simplification will be obtained. It la possi- 
ble that propagation angles may be easier to mea- 
sure than cut-off rotioB, although either one would 
be sufficient, ThiB problem Is currently being 
atudled, 
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INCIDENCE ANGLE. <P r . deg 

Figure 2. - Correlation of optimum resistance with wave 
incidence angle on the wall, tj • 15. M Q • -0.4, e • 
0.05. 




INCIDENCE ANGLE, <p r deg 

Figure 3. - Correlation of optimum reactance with wave in 
cidence angle on the wall, n ■ 15, M 0 • -0.4, t ■ 0.05. 







Figure 4. - Effect of boundary layer thickness upon optimum 
resistance-incidence angle correlation, r? ■ 15. M 0 ■ -0.04. 
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Fiqure 5. - Effect of boundary layer thickness upon optimum 
reoctance-inciJence angle correlation, 7j • 15. M 0 • -0.4. 
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Figure 6. - Correlation of optimum reactance with resultant 
incidence anole. £ • 0.005. tj ■ 15. M 0 • -0.4. 
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Figure 7. - Sketch of wave front and propagation angles in a rec- 
tangular duct. 



MAXIMUM ATTENUATION 

Figure 8. - Comparison of approximate and 
exact calculated attenuations over the 
complete attenuation range, rectangular 
duct, second symmetric mode, M ■ 0. 

V ' 5. 
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RATIO OF EXACT ATTENUATION TO EXACT 
MAXIMUM ATTENUATION 

Figure 9. - Comparison of approximate and 
exact calculated attenuations over the 
complete attenuation range, rectangular 
duct, second symmetric mode. M ■ -0.4. 
V -5. 
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1C) EIGENVALUE MAGNITUDE. 

Figure 10. - Propagation calculations made 
over entire range of propagation coeffi- 
cient. i. rectangular duct. M • 0. o • 

0. 02. t? ■ 5. 
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RATIO OF EXACT ATTENUATION TO EXACT 
MAXIMUM ATTENUATION 

W Figure 11. * Comparison of approximate and 

exact calculated atienuations over the com 
plete attenuation range, cylindrical duct, 
seven lobed spinning mode, second radial, 
M • 0, n • ?0. 


